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Abstract. Here we present the results of applying the generalized Riemann £- 
function regularization method to the gravitational radiation reaction problem. 
We analyze in detail the headon collision of two nonspinning black holes with 
extreme mass ratio. The resulting reaction force on the smaller hole is repulsive. 
We discuss the possible extensions of these method to generic orbits and spinning 
£\J , black holes. The determination of corrected trajectories allows to add second 

■ perturbative corrections with the consequent increase in the accuracy of computed 

' waveforms. 

o 
o 

PACS numbers: 04.30.-w, 04.25. Nx, 04.25.-g, 04.70. Bw 

o 
o 

""^ ' There has been increasing astronomical evidence during the last thirty years in 

q-i favor of the hypothesis that supermassive black holes, lying in the center of galaxies 

are a common phenomenon in our Universe. They are thought to be the main engines 
(3JT). of quasars, X-ray sources in the core of galaxy super-clusters, and gamma ray bursters. 

There is even compelling evidence that our own galaxy shelter a black hole of three 
million solar masses in its center Although compelling, this is still indirect evidence 
of the existence of black holes. It appears that their characteristic gravitational wave 
signature provides the best way to determine its definitive existence. The LISA project 
j| offers an excellent opportunity to match theory and observation in order to reach 
such goal. From the characteristics of this space detector, the main astrophysical 
sources expected to be "seen" are precisely supermassive black holes accreting stars 
close to the galactic center. This scenario offers the possibility to use perturbation 
theory around a single black hole as the appropriate theoretical model to compute 
gravitational radiation from such extreme mass ratio binary systems. 

Perturbation theory reached maturity already in the early 1970's|}| ||. Linearized 
Einstein equations have been successfully combined into wave equations for the two 
degrees of freedom of the gravitational field propagating in the curved background of 
a massive black hole. The source terms depend on the trajectory of a test particle 
moving along a geodesic of the background metric. There remained one important 
problem to be solved: The changes induced in the trajectory of the particle due to the 
radiation generated by the orbital motion of the same particle around the massive hole. 
Progress along this line always stumbled with the infinities appearing in the equation 
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of motion of the particle when perturbations are evaluated, precisely at the location of 
the particle. A consistent regularization method was lacking until recently, when Mino, 
Sasaki and Tanaka|5| provided a renormalized equation of motion using two different 
methods: i) The Hadamard regularization and ii) An asymptotic matching of the 
near and far metric expansions. Quinn and Wald |6) assume an axiomatic approach to 
obtain independently the same equation of motion. While the above work provide well 
founded formal expressions, they appear yet impractical when dealing with concrete 
computations since they rely on the use of the harmonic gauge. Here we present 
an alternative, practical regularization method based on the analytically continued 
Riemann function. 

A test particle with rest mass too is here represented by the energy-momentum 
tensor 

TJI1TJV 

T^ = m 01J ^S[r-r p (t)}S^^ P } • (1) 

Since the only "forces" acting on the particle are gravitational, the equation of 
motion of a test particle is given by the geodesic equation 

It is precisely when one wants to evaluate the Christoffel symbols at the location of the 
particle that divergences appear. They originate in the Dirac's delta representation 
we give to the source term ([!]) of Einstein equations. 

In Refs. §, § 1 we studied the headon collision of two Schwarzschild (nonspining) 
black holes with masses M and m respectively. We treat this problem as perturbations 
of the Schwarzschild background. The only expansion parameter being to/M, assumed 
to be much less than one. The symmetry of the background allows us to decompose 
metric perturbations into spherical harmonics Yi m (9, (f>) and the gravitational degrees 
of freedom satisfy the wave equation 

d 2 ijj e d 2 ip e T/ 1 \ i c ( n 
" -Qp + q^2 ~ = S ^ *) ' ( 3 ) 

where r* = r + 2M ln(r/2M — 1), Vt is the Zcrilli potential, and Se(r, t), is the source 
term generated by the small hole J?], ||. 

From the waveform ipg(r,t) one can reconstruct all perturbed metric coefficients 
in the Regge - Wheeler gauge. One can then prove that all metric coefficients are 
C° at r p for every I. The divergencies reappear when trying to add up all multipole 
contributions To see this more explicitly, let us consider the r and t components 
(the only nontrivial for the headon case) of the geodesic equation (||), and combine 
them into a single equation for r(t) 

§ = ^ (£V+(n- - r; r ) - 2Ti) ( ^Vn,(4) 



dt 2 rr \dtJ y tr rr '\dt) yu tr '\dt) 

Linearization of this equation and subtraction of the zeroth order geodesic gives the 
deviation of the trajectory Ar p (i) from the zeroth order one 



r p (t) = d tZp = -(1 - 2M/z p y 2 ^YM/T° ' (5) 
directly in terms of Schwarzschild coordinates, ready for further applications 

Af p = A Ar p + B Ar p + C (6) 
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It is precisely Cg the piece that being finite for each multipole diverges when summed 
over £ (here we consider the averaged value of Ct over both radial sides of r p ). The 
direct numerical integration of the wave equation shows that we can split Ci as 



C = 



OO 

£ {(2^ + 3/2 - 0)-' 3+1/2 C oo + Cr n ) , 



(8) 



where the introduction of /? is motivated by the D-dimensional extension of the initial 
value problem for time symmetric conformally flat data. In this case —j3 = D — 9/2. 
We neatly determined numerically that (3—1/2 along the trajectory of the particle. 

This form of the coefficient C, determining the corrected trajectory of the particle, 
shows clearly that the I— independent term, summed over £, diverges. We also 
determined numerically that C£ en behaves like £~ 2 for large £, thus giving raise to a 
finite contribution to C. The key observation here is that we can bring the divergent 
sum to the form of the Riemann (— function [|ll], ((a, b) = YlT=o^- + b)~ a , 



C = 2- /3+1 / 2 C 0O C(/3 - 1/2, 1/2) + 



CI 



(9) 



Since the analytically continued ^-function gives £(0, 1/2) = 0, in order to regularize 
C, we must just subtract to each multipole their £ — > oo piece. When we do that 
[fiol we have found that the effect of C ren on the trajectory is to generate a repulsive 
contribution that has a maximum near the pick of the Zerilli potential at r = 3.1M. 

Let us label functions with the ± superindex to refer to their values in the region 
r > r p and r < r p near the location of the particle r p respectively. One can see that 
both, initial data and Zerilli equation imply the following formal symmetry 

s/Z<tpf(-L) = Vlijf(L) (10) 

where L — £ + 1/2. The same symmetry holds for the t and r derivatives, hence the 
reaction 'force' will have the form ^ r± (L) = p(L 2 )±L q(L 2 ), and the large L expansion 
of its average 



A,: 



£ = 



e=o n=0 



71=0 



A, 



1=0 



2n. 



£ 

n=0 



A n C(2n, 1/2) 
A C(0, 1/2) + A 2 C(2, 1/2) + Ai C(4, 1/2) 



(11) 
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where we observe that £(0,1/2) = gives the regularization, £(2,1/2) = i w 2 gives 
the leading reaction 'force' and £(4, 1/2) = 1 7r 4 gives an estimate of the ^ error 3 of the 
leading term of the order of 7 — 8%. 

Barack, Ori and Burko Jl3|, [l4| have developed independently a similar 
approach (in the sense that they also decompose into multipoles), but that uses 
the regularization method of Ref. |Bj and implements it to scalar radiation in the 
Schwarzschild background. Their results are completely compatible with those of the 
£— function regularization in the sense that in the end it is exactly only the I — > oo 
that has to be subtracted to the originally divergent force. 

With Barack Jljj we have been able to confirm analytically the numerical behavior 
discussed above. Based on a local analysis of the two-point Green's function we make 
a 1/L (with L = I + 1/2) expansion to determine the behavior of the waveform ipg 
and its derivatives for large £, for instance 



ipf(r p ) =4m V2^v r L{L- 3 ±2EL~' i + O(L~ 5 )} 
d r tp e (r p ) 



± = W5jF^f a _ 3 

" . 2.1/ I 2 



dti>t(r p ) = - (1 - 2M/r p ) (^j d r ^(r p ) (12) 

for a particle released from rest at ro, E = y/l — 2M/r . ± stands for the side 
derivatives about r p (waveforms are discontinuous). Hence we can compute 

Cf = a ± L + b + c ± L- 1 +0(L- 2 ), (13) 

where a^, b, and are Z-independent coefficients. We find that the average value of 
a and c vanish at the location of the particle, while the average 

<C(l^oo)>=-^-^(l-2M/r p f/ 2 E- 2 (14) 

This method can be used to assist the numerical computation for large I carrying 
out the expansion to order 1/L 2 in C. Thus, making necessary the numerical 
integration for only the first few lower multipoles. This might be crucial when dealing 
with generic orbits around Schwarzschild black holes. The analytic expansions would 
allow a clean application of the £— regularization scheme and sum over large £ of the 
regularized reaction " force" . 

In the orbital case it will also be important to test our final results with those of 
the energy-balance approximation Jll], [ll| |l7j| This method provides accurate results 
for circular orbits on the Kerr background |lq| . Our method, in principle, could also 
be implemented to a particle orbiting a Kerr hole since the Teukolsky equation, that 
describes the perturbations around rotating holes, can be decomposed into multipoles 
in the frequency domain (Laplace decomposition of the time-dependence [0). The 
details of this implementation remain to be explicitly worked out. 

A further application of the above results one can foresee is the extension of the 
analysis to second perturbative order JTgj] . In this way one can increase the accuracy 
with which waveforms are computed as well as reaching not so small mass ratios m/M. 
At this stage, the resulting waveforms will not only be relevant for LISA, but also for 
ground based gravitational wave detectors sensible to galactic binary black holes of 
comparable masses and black hole - neutron star systems. 
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I would like to end this contribution with an optimistic note from the 
theoretical side and predict that the subsequent progress in our understanding of 
radiation reaction will soon bring good news for modeling gravitational emission from 
astrophysical sources. 

Acknowledgments: Thanks to L. Barack for many useful remarks. 
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Appendix A. Connection coefficients 



Here we make explicit the details that enter into the computation of Ci given in Eq. 
(0). The connection coefficients involved in the first order geodesic, Eq. (|J) in terms 
of the first order perturbations as defined in the Regge - Wheeler (RW) gauge || are 
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where the overdot and prime mean derivative with respect to the Schwarzschild time 
t and radial coordinate r respectively. In addition, we have from Einstein equations 
in the RW gauge (see Ref. jj], Eq. (C7g)), that for a radial infall 

H = H 2 , £>2. (A.2) 
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This allows us to write everything in terms of Hf and H% only (and its derivatives) . 
K l does not appears in the headon geodesies. 

Proven the continuity of the metic coefficients at r = r p our way is clear to 
compute the trajectory of the particle to first perturbative order. Our last step is 
to compute these metric derivatives in terms of tp and ip. Since we know that the 
metric coefficients are continuous, their derivatives will have a jump, but not delta 
(or derivatives of it) terms. Then, in order to simplify expressions we will not write 
the source terms and use an overbar [ip) on radial derivatives of ip to indicate that we 
have already subtracted its singular behavior at r = r p and taken the average value. 

The r - derivative of Hi from Ref. (HJ is 

n » ,n T r Ar 2 - XMr + 3M 2 ~ 
d r H 2 = (r - 2M) d?il> ■ ■ 



1 



-A 2 (A + 1) r 3 — A (2A - 3) Mr 2 



(Ar + 3M)r 
- 15AAf 2 r - 



d 2 r V 



18M £ 



r 2 (\r + Miy 



d r ip 



(A.3) 

A 3 (A + 1) r 4 + 3A 2 (A - 1) Mr 3 + 27A 2 Af 2 r 2 + 63AA/ 3 r + 54M 4 



r 3 (Ar + 3M )" 



4> 



Likewise, upon r 
L0| we obtain 



derivation of expression the expression for Hi, given in Ref. 
(3 - 5A) Mr 



d r Hi = rd 2 ip + 



2Xr 2 



[A (A 



9AP 

(r - 2M) (Ar + 3M) 
- 3) r 2 - 6AAfr + 3 (4A -f 



<9 r ^ 
3) M 2 ] M 



(r - 2M) (Ar + 3M) 2 
The time derivative of Hi can be also obtained 



(A.4) 



9M 3 + 9AAf 2 r + 3A 2 Mr 2 + A 2 (A + l)r 3 
r 2 (Ar + 3M) 2 



ip (A.5) 



and finally we find dtHi from Zerilli'sQ E q. (C 7e), where K' can be found from 
Hi, given in Ref. |10| , and #2 is given by Eq. (A.5), 

(r-2Af) 2 , - , , (l5M 2 + 3AMr + Ar 2 ) „_ 

dtHi = i ^9 3 V + r - 2M ^ 2 , . , ^ 2 </> 

r H (3 Ai + A r) 

- 2 (2 A 2 A/r 3 + 18 A A/ 2 r 2 + 27 rM 3 - 9 A 2 Af V - 42 A Af 3 r 

- 63 M 4 + A 3 r 4 + A 2 r 4 - 2 MA 3 r 3 ) / (r 3 (A r + 3 Af) 2 ) d r i> 

+ 2 (-189 M 5 + 81 Af 4 r - 216 A A/ 4 r + 90 A A/ 3 r 2 - 90 A 2 AT 3 r 2 

+ 33A 2 Af 2 r 3 - 18A 3 Af 2 r 3 - 3A 3 (A - l)Afr 4 + A 3 r 5 + A 4 r 5 ) ^ 

/ (r 4 (3M + Ar) 3 ) (A.6) 

These expression complete the equations needed to (numerically) integrate the 
geodesic equation (0) that gives us the particle's trajectory to first perturbative 
order. Note that the numerical implementation should be able to handle third order 
derivatives of tpf. 



